Optical conductivity in non-equilibrium d-wscve superconductors 
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We consider the optical conductivity of a d-wave BCS superconductor in the presence of a non- 
equihbrium distribution of excess quasiparticles. Two different simplified models used in the past 
for the s-wave case are considered and results compared. In the r*-model of Parker the excess 
quasiparticles are assumed to be in a thermal distribution at some temperature T* larger than the 
equilibrium sample temperature. In the /i*-model of Owen and Scalapino a chemical potential is 
introduced to accommodate the excess quasiparticles. Some of the results obtained are specific to 
the model, most are qualitatively similar in both. 
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I. INTRODUCTION 

The study of non-equilibrium superconductivity in 
conventional s-wave superconductors has a long history. 
Several comprehensive reviews exist iSi^ Much less has 
been done for unconventional order parameters includ- 
ing the high Tc oxides. The cuprates are particularly 
interesting since there is now a consensus that its energy 
gap has (i-wave^i^i^ rather than s-wave symmetry. This 
difference is significant particularly since it implies the 
existence of gap nodes on the Fermi surface and hence 
the existence of quasiparticle states of low energies. 

One might expect that quasiparticle excitations cre- 
ated at rather high energies through application of a 
laser pulse to the equilibrium sample or by some other 
means, would rapidly relax towards the lowest energy 
states available and accumulate in the nodal regions of 
the gap on the Fermi surface where the density of excited 
states goes linearly in energy. The fast relaxation pro- 
cesses can proceed through the electron-electron interac- 
tion which is strong in the cuprates and also through cou- 
pling to the phononsiS The phonons are a separate system 
of excitations which can share in the excess energy pro- 
vided by the excitation source. This is to be contrasted 
to what happens in formulations of highly correlated sys- 
tems in terms of boson exchange processesj^iiSiiiii^iiii^ 
such as the exchange of spin fluctuationsi^ in the case 
of the Nearly Antiferromagnetic Fermi Liquid (NAFFL) 
model. While such exchange processes can equilibrate 
the energy, they do not take it out of the electronic sys- 
tem. Another feature of d- as opposed to s-wave non- 
equilibrium superconductivity is that, as has been argued 
by Feenstra et al.^ the final stage of equilibration which 
proceeds through recombination of nodal quasiparticles 
into Cooper pairs, may take a much larger time than in s- 
wave because of restrictions on the kinematics. New fast 
laser pulse techniques^iiSjiS have helped greatly in the 
study of non-equilibrium state in cuprates. These tech- 
niques have recently been used in leadi^ to measure the 
conductivity (j{T, w) under non-equilibrium conditions as 



a function of frequency in the entire infrared range. Simi- 
lar experiments in the cuprates should be possible and we 
hope this paper will stimulate such work. Here we pro- 
vide a calculation of non-equilibrium effects just before 
recombination on the optical conductivity of a d-wave su- 
perconductor. We also consider explicitely the microwave 
region where analytic results can be obtained. 

Very recently Nicol and Carbotte^S studied the ef- 
fect of a non-equilibrium distribution of quasiparticles 
on the gap in a BCS d-wave superconductor. In this 
work two well established models coming from the liter- 
ature on conventional non-equilibrium superconductors 
are employed. These are the /.t*-model by Owen and 
Scalapinc-^^ and the r*-model of Parker i22i2^ These mod- 
els, while very simple, have been useful in conceptualizing 
what can happen, in planning, and interpreting experi- 
ments. In some sense they represent two extreme limits 
between which the true distribution might be expected 
to fall. The /^*-model assumes that the excess quasiparti- 
cles rapidly accumulate in the lowest quasiparticle energy 
states available which is around the gap A in s-wave but 
is around zero in d-wave where the density of states shows 
a linear in energy dependence. A finite chemical potential 
/i* (for the quasiparticles) is used to describe the occu- 
pation. The other extreme is to assume that the excess 
quasiparticles equilibrate at some new temperature T* 
larger than the sample temperature T. In both cases the 
system is not in thermal equilibrium because the mea- 
surement is to be made before the final recombination 
process into Cooper pairs proceed. This is assumed to 
be the largest time scale in the relaxation process. 

There are significant differences in behavior for d-wave 
as compared to s-wave. For example in terms of the 
excess quasiparticle density n (to be defined more pre- 
cisely later) the reduction in the gap (5A/A is given by 
—2n for both, T* and ji* models in s-wave, but for d- 
wave it goes like — (4-\/2/3)Vn^ for the /i*-model and 
-(32/7r3)^(3n)3 in the T* model. These fundamentally 
different dependences on n, linear in s-wave, and to a 3/2 
power in d-wave can be simply understood. The presence 



of excess quasiparticles blocks states around the Fermi 
surface which would otherwise be available to form the 
condensate in a variational sense. In d-wave, the exci- 
tations are assumed to gather around the nodal regions. 
In these regions the gap is small and contributes little 
to the condensation energy so that blocking such states 
is less effective than in the s-wave case where the gap is 
everywhere finite and equal to A. 

In their paper Nicol and Carbotte^ii consider ex- 
plicitly the case of tunneling in a normal-insulator- 
superconductor (N-I-S) junction with a non-equilibrium 
distribution on the superconducting side. They also ad- 
dress the case of pump probe experiments which claim to 
measure the temperature dependence of the excess quasi- 
particles n{T) before final recombination into Cooper 
pairs. In the present paper we extend the work to the 
optical conductivity. Such calculations exist in the s- 
wave case^"* based on an appropriate generalization of the 
dirty limit Mattis-Bardeen formula^^ for the optical con- 
ductivity. Here we consider the d-wave caseSS and treat 
explicitly both, the T* and the fj,* model and compare 
results obtained for the same value of excess quasiparti- 
cles n. Since in many aspects both models give the same 
qualitative physics, although there are some important 
differences which we will note, we can expect that more 
realistic distributions would not be so different. 

In Section II we introduce some elements of the the- 
ory of the optical conductivity cr{uj,T) as a function of 
frequency uj at temperature T suitably generalized for 
the non-equilibrium case. In the T*-model no formal 
changes to the formulas are needed, only the tempera- 
ture is to be interpreted differently. In the //*-model all 
necessary changes to the mathematical formulas can be 
incorporated through a change of the thermal factors. In 
Section III we present exact numerical results for the fre- 
quency dependence of the optical conductivity (real and 
imaginary part) as well as for the reflectance and the 
normalized reflectance difference due to an excess dis- 
tribution (n) of non-equilibrium quasiparticles which is 
the quantity that is directly measurable. Results for the 
T* and /i*-model for the same n are compared. Section 
IV contains a number of analytic results for a regime in 
which the temperature is larger than the impurity scat- 
tering rate i.e. the weak scattering limit. Comparison 
with exact numerical results is also given. Dependences 
on the excess quasiparticle density n are made explicit 
and helps physical understanding. A short conclusion is 
found in Section V. An Appendix contains some mathe- 
matical details. 



the introduction of a chemical potential fj,* and in the 
other by simply changing the temperature from T, the 
temperature of the sample before application of a laser 
pulse, to T* the non-equilibrium quasiparticle tempera- 
ture before final recombination. For simplicity, we will 
assume here that the initial sample temperature T is suf- 
ficiently low that it can be considered to be zero. In the 
actual numerical work it is taken to be 3 K for conve- 
nience with a gap of 24meV. For a BCS d-wave super- 
conductor this choice corresponds to a Tc of 130 K. If we 
further assume that the density of excess quasiparticles n 
is small, we can obtain analytic results for the change in 
gap and for the corresponding value of /i* and/or T* . The 
necessary results are to be found in the paper by Nicol 
and Carbotte.2^ In the r*-model n = {n'^ /12){T* /Aof 
and {SA/Aq) = -{32/TT^)^/(3nj^, while in the /x*-model 
fi* /Ao = V2^ and SA/Aq = -{AV2/3)Vri^. [From 
now on we will use A or A{n) for the gap with a fi- 
nite n and Aq for its equilibrium value when n = 0, 
i.e. Aq — A{n — 0).] These relationships can be used 
directly in our conductivity calculations. Some expla- 
nation is needed. The above results were derived for a 
d-wave BCS superconductor in the pure limit. To get 
meaningful conductivity results it is necessary to for- 
mally include some scattering mechanism. Here we will 
simply treat elastic impurity scattering. When this is 
done the impurities also formally enter the gap equation 
and in principle one should include these modifications 
along with those brought about by the non-equilibrium 
distribution of quasiparticles in the formula for the con- 
ductivity. Here, for simplicity, we will not do this. The 
impurities do reduce the size of the d-wave gap ampli- 
tude but we will assume that this is already included in 
our initial choice of this parameter and simply treat the 
changes brought about by n based on the pure case. As- 
suming that the impurity concentration is very small as 
we wish to do here, this is sufficient. There will also be 
a corresponding effect of impurities on the chemical po- 
tential fi* but this is also neglected in the limit of small 
impurity scattering. With this explanation we can pro- 
ceed with the calculation of the optical conductivity. 

It is conventional, and we will follow this definition 
here, to measure the excess quasiparticle density n in 
units of 4N{0)Ao where the factor 4 is from spin de- 
generacy and restricting the defining integral to positive 
energies only. iV(0) is the normal state electronic density 
of states at the Fermi surface. In these units 
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II. THEORY 

The non-equilibrium distribution in either T* or fi*- 
model does not change the symmetry of the gap. It does 
reduce its magnitude. In addition, the thermal factors 
which enter directly the infrared conductivity formula are 
also changed. These changes are modeled in one case by 



for the /^*-model and 
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in the T*-model. Here /T(e) is the usual Fermi Dirac dis- 
tribution /T(e) = [l + exp(/3e)]~^ with fi = l/fc^T where 



ks is the Boltzmann factor which we set to one from here 
on. In Eqs. ^ the brackets {■ • •)e imply an angular av- 
erage over the polar angle 9 which gives the position on 
the Fermi surface in the two dimensional Cu02 Brillouin 
zone. In terms of this angle the gap A{9) = Acos(26'). 
The quasiparticle energies -Ek — -^/e^ + A^(r, 9) where 
we make explicit the temperature dependence of the gap 



which enters the T*-model. Here k is momentum and 
eic is the normal state band energy. Consistent with our 
treatment of the gap equation the clean limit is used in 
writing down Eqs. ^. 

The formula for the optical conductivity in the equi- 
librium case takes on the formiSSiSLSSiSS 



J 



2 / °° 



where Vlp is the plasma frequency and the function J{lu, i>) is 
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[1 - N{lu)N{uj + v)- P{uj)P{uj + v)] 
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with i?i (w), N*{uj), and P*{uj) the complex conjugate of 
Ei{lo), N{u!), and P(w), respectively. In Eq. (jSJl 



Ei{uj) = \Jco^{ijJ + iO+) - A2(tj + iO+), 
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In these equations uj{uj + iO^) is the renormalized fre- 
quency which includes the impurities and in the pairing 
energy A{uj + iO~^) the angular dependence on 9 has been 
suppressed. a'{T,v) is the optical conductivity as it is 
calculated within our numerical programs. It is given in 
units of meV~^. If Q,p is given in meV then a(T,v) is 
given in units of meV which can easily be transformed 
to SI units using the relation that one Vl~^ va~^ corre- 
sponds to 5.916 X 10^3 meV. Finally, in BCS A(tj-hiO+) 
is independent of uo. 

To treat a non-equilibrium system a chemical poten- 
tial ^* for the excess quasiparticles is introduced in the 
grand canonical ensemble average and this means that 
the thermal factor in Eq. Q needs to be replaced by: 
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while in the T*-model there are no modifications except 
to change the temperature T to T*. In Fig. ^ we com- 
pare the two factors: F{uj) and ta,nh{uj /2kBT*) with /i* 



related to T* by T* ~ 1.56 ^*. We see that both are anti- 
symmetric with respect to uj = and that F{lu) depletes 
the region around uj = more than does the hyperbolic 
tangent. These differences will reflex themselves in the 
conductivity. We begin by giving results of our numerical 
evaluation of Eq. (|2Jl. Before proceeding, it is necessary 
to specify how impurities are to be included in the theory. 
The impurities renormalize the frequencies to uj^lu + iO^) 
from uj according to the equation22i2i 



u{uj + iO+) = w + inV^ 
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where we have introduced the impurity scattering in a 
T-matrix formalism, c = 1/ {2TrN{0)Vimp)- Here Vimp is 
the strength of the impurity potential and A''(0) the elec- 
tronic density of states at the Fermi surface in the normal 
state which is taken to be constant, i.e. energy indepen- 
dent in the range of energies relevant to superconductiv- 
ity. In Eq. Q) the real part of {N{oJ))g is the quasiparticle 
density of states. The parameter TrF^ is related to the 
impurity concentration nj by ttF^ = n//[iV(0)7r]. When 
Vimp is very small c — > 00, we recover the Born limit 
in which case we denote F^/c^ as t"*", and for Vimp very 
large c — > which is called the unitary limit. A realistic 
case would be for finite c. Fits to experimental data have 
given values of the order of 0.2jS, From the definition of 
c and taking 7V(0) to be 1/W for c = 0.2, Vi„,p ~ 1.3 W, 
where W is the band width. 
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FIG. 1: The thermal factor tanh{uj/2kBT*) which en- 
ters the optical conductivity (Eq. ||5J) in the T'-model 
(solid curve) compared with the equivalent function F{uj) = 
[—6{u} + n*) — 9{ui — n*) + 1] which enters the /x* model 
(dashed curve). As shown in the text for a common value 
of excess quasiparticles n, 2T* = (2\/6/7i")A'* — 1.56/^*. This 
relation is built into the figure. Here T* = 53 K and ^* — 5.88 
have been chosen. 



III. NUMERICAL RESULTS 

In the top frame of Fig. 13 we show results for the 
real part of the optical conductivity (j[(uj) as a func- 
tion of frequency w in computer units. To get the actual 
conductivity one must multiply the computer results by 
flp/An with ftp the plasma frequency. The parameters 
used for the run are Aq = 24 meV, T = 3 K which is low 
enough that it is representative of zero temperature and 
t+ — 0.1 meV. The impurity scattering is taken in the 
Born limit (c — > oo). The solid gray curve labeled simply 
by Aq on the figure is the conductivity in the equilib- 
rium limit, i.e. n = 0, and is included for comparison. 
The other four curves are with a finite non-equilibrium 
excess quasiparticle density of n = 0.03. This value of 
n is large compared to a value that one may realistically 
have in an experiment. The results obtained, however, 
are illustrative of what one might expect even for smaller 
n. For the /j,*-model /i* ~ 0.25Ao and the change in 
the gap (5A ~ — O.OIAq which is very small. In the T*- 
model the corresponding number is T* ~ 0.19Ao and 
SA = — 0.028Ao, again small but three times bigger than 
in the /i* model. In the top frame of Fig.|21the solid black 
and dashed curves are for T* = 53 K with and without 
inclusion of the small gap change. We see that the direct 
thermal effects of T* are much more important than any 
slight change in the gap. This also holds for the /x*-model 
with dotted (no gap change) and dash-dotted (with gap 
change) curves. Comparison with the n — case shows 
that the real part of the conductivity is strongly affected 
by the inclusion of a chemical potential /i* = 5.88 meV 
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FIG. 2; Top frame: the real part of the optical conductivity 
(^[{uj) as a function of frequency uj. The solid gray is for the 
unirradiated equilibrium sample with temperature T = 3 K, 
gap Ao = 24meV, and an impurity content described by 
t^ — 0.1 me V. The solid (dashed) curve includes an excess 
quasiparticle density n = 0.03 in the r*-model (T* = 53 K) 
with the gap unaltered (altered to 0.972Ao). The dotted 
(dash-dotted) curve is for the /^*-model {fi* = 5.88 meV) with 
the gap unaltered (altered to 0.99Ao). The bottom frame 
is the same as the top frame except that now the impurity 
content is greatly reduced to i+ — 0.01 meV, with a single 
T* = 58K (solid line) and a, fi* = 5.88 meV (dashed line). 



and that there is significant difference between the pre- 
dictions of the T* and /i* models. In particular, for the 
/i*-model there is a distinct structure predicted to occur 
at the value of the gap Aq minus the chemical poten- 
tial fj,* . The characteristic of these structures is more 
easily seen in the bottom frame of Fig. |2] where addi- 
tional results for (j[{i-j) vs lu are given for the same gap 
value of 24 meV but now the impurity content is reduced 
by an order of magnitude to f+ = 0.01 meV. The gray 
solid curve which applies to the equilibrium case, shows 
a slight change in slope at w = Aq = 24meV as well as 
at w = 2Ao = 48 meV. These structures are characteris- 
tic of the logarithmic van Hove singularity in the d-wave 
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FIG. 3: Top frame: same as in the top frame of Fig.|21except 
that now uia'2{u!) is plotted, i.e. uj times the imaginary part of 
the conductivity which is related to the inverse square of the 
penetration depth at a; ^ 0. Middle frame: the reflectance 
as a function of to. Bottom frame: the normalized change in 
the reflectance 5R{ui)/ R{<jj) vs uj. To get these we have used 
fip = 2eV and Soo = 1. 



quasiparticle equilibrium density of states. In compari- 
son to the gray solid curve, the dashed curve which has 
/i* — 5.88 meV shows an additional pronounced structure 
at ~ 18 meV which is the value of Aq — n* in this case. 
When we also include a small shift in the gap we do see 
it as a small shift in the 2Ao structure which corresponds 
to the smallest optical frequency which can connect both 
logarithmic van Hove singularities in initial and final elec- 
tron state. By contrast, the structure at Aq — ^* corre- 
sponds to a transition between the top of the occupied 
region of the density of states as the initial state with 
the van Hove singularity as the final state. The remain- 
ing curve in the bottom frame of Fig. |21 is for the T* 
model (solid line). Now the region below 2Ao is quite 
smooth and shows no discernible structure in contrast to 
the /x*-model. This is expected since in the r*-model the 
non-equilibrium distribution is assumed to be thermal. 

Fig. O shows additional results for the case f^ — 
0.1 meV. In the top frame we show Lija2{u>) (w multi- 
plied by the imaginary part of the conductivity) also 
denoted by a frequency dependent inverse penetration 



depth A^ i'-^^T) in the literature. The middle frame 
shows the reflectance while the bottom frame shows 
the normalized difference in reflectance between non- 
equilibrium and equilibrium case, i.e. 5R{uj)/ R{lo). This 
is often the quantity that is measured directly in non- 
equilibrium pump probe experiments and is included here 
for convenience. Considering the top frame first, we see 
significant differences in the value oi Loa2{Lo) for frequen- 
cies Ld below the gap Aq — 24 meV. In particular, a peak 
which does not exist in the solid gray curve for the equi- 
librium case develops in both, the T* and the /i*-model. 
The dotted and dash-dotted curves apply to the /i*-model 
without and with change in the gap included while the 
solid and dashed curves are for the T*-model. The the- 
ory predicts a sharper peak which forms at smaller en- 
ergy for /i* as compared with the r*-model. The center 
frame of Fig. O gives the reflectance R{lo) vs lu for five 
cases as in the other frames. As in the top frame, the 
largest differences in this set of curves occur at frequen- 
cies below the gap to < 24meV. In that case the T*- 
model (solid black and dashed curves without and with 
the non-equilibrium change in the gap included) predicts 
a very rapid drop in R{uj) at the lowest values of co and 
then a plateau before a second rapid drop which sets in 
around the gap. The curve for the /i*-model, dotted and 
dash-dotted without and with gap change, is in compar- 
ison more structured and in fact shows a small peak in 
the region 10 to 15meV. The bottom frame of Fig. ^ 
serves to emphasize these differences even more. What 
is plotted is the difference between non-equilibrium and 
equilibrium value of R{uj) normalized to this equilibrium 
value. The quantity —5R{uj)/ R{uj) shows a very rapid 
rise out of cj = which is larger for the T* model than 
for the /i*-model (by roughly a factor of four). After a 
maximum is reached, a minimum is seen only in the /i* 
case which shows a second maximum within the gap re- 
gion (corresponding to the structures in the conductivity 
in the top frame of Fig. [3 at Aq — /x*) before the change 
in reflectance becomes quite small. In both models the 
region between Aq and 2Ao is also affected and is quite 
sensitive to changes in the gap value. 

We have carried out additional calculations for val- 
ues of c in the impurity potential away from the Born 
limit. Here we report only one case. In Fig. ^ we show 
results for Aq — 24meV, as in a previous figure with 
T = 3K, r+ = 0.01 meV and c = 0.2 which is near the 
unitary limit. What is shown in the top frame is the 
reflectance as a function of frequency. The gray curve 
is the equivalent equilibrium case included for compar- 
ison. The dashed curve is based on the r*-model and 
the dash-dotted curve on the /i* -model. The middle 
frame gives the corresponding results for the real part 
of the optical conductivity. Comparison of these results 
with the equivalent results shown in the bottom frame 
of Fig.|21reveales that the Born limit is a better case in 
which to investigate non-equilibrium effects in the sense 
that, for the same non-equilibrium quasiparticle density 
n the effects are less pronounced for the c = 0.2 case. 



IV. SIMPLE ANALYTIC RESULTS 
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FIG. 4: Reflectance R{uj) (top frame), real part of the op- 
tical conductivity (7'i{lu) (middle frame), and normalized re- 
flectance difference 5R{lo)/ R{lo) (bottom frame) vs uj for a 
case with Ao = 24 meV, T = 3 K, r+ = 0.01 meV, and c = 0.2 
near the unitary limit. The gray curve in each of the top two 
frames is the equilibrium case shown for comparison. The 
dashed curve is for the r*-model and the dashed-dotted curve 
for the /x*-model. 



In contrast to the Born limit, the <j'i{uj) vs lj (solid gray 
curve) is now much closer to the non-equilibrium results 
(dashed or dash-dotted curves). The normalized differ- 
ence 5R{uj)/ R{u!) vs w is shown in the bottom frame 
of Fig. 21 The differences between the r*-model (dashed 
curve) and the /i*-model (dash-dotted curve) are less pro- 
nounced than was found in the bottom frame of Fig.|21 In 
our previous analysis of the optical conductivity of vari- 
ous cuprates we have found, in several cases, the need to 
include in our calculations some impurity scattering in 
the unitary limitti^i^S On the other hand, in some of the 
very pure samples grown in BaZnOs crucibles the mean 
free path at low temperatures is found to be of the order 
of a micron and a value of c 7^ is needed J^ 



It is possible and instructive to derive some simple an- 
alytical results for the non-equilibrium conductivity cal- 
culated numerically in the previous section. This can 
provide insight into the physics involved as it makes ex- 
plicit the dependence on excess quasiparticle density n of 
various quantities. In terms of the London penetration 
depth A^^(O) the penetration depth in the pure limit at 
finite temperature and finite /i* is given by 
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(5) 
where N{E) is the quasiparticle density of states. For 
small /I* at zero temperature the thermal factors in 
Eq. ini become a (S-function; we can also use the nodal 
approximation for N{E) — E/Aq and get immediately 



K' 



{0,fi*) = Xf{0) 



1-f 

Ao 



- X-2 
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The London penetration depth at zero temperature is in- 
creased by the presence of a non-equilibrium distribution 
and the reduction in superfiuid density follows a square 
root of n law. This is also seen in the results presented 
in the top frame of Fig. 3. A similar result also holds for 
the T* model. From Eq. © with /i* = and T replaced 
by T* we get 



Xf{T*,in = AZ'(O) 



Ai'(O) 




(7) 



The intermediate formula in 10 is just the well known 
linear in T law for a pure d-wave superconductor. The 
coefficient of the ^/n term in the last expression in Q is 
1.53 while in the /i*-model the coefficient was 1.4. For the 
same non-equilibrium density n the superfiuid density is 
less reduced from its equilibrium value in the /i*-model. 
This mirrors but is the opposite of what was found for 
the superconducting gap. In both models the reduction 
follows a VTV^ law with coefficients 1.9 in the ^*-model 
as compared with 5.4 in the T*-model. The physics of 
the reduction in this case is that the states around the 
Fermi surface occupied by the excess quasiparticles more 
effectively block the formation of the condensate in the 
/i*-model than they do in the r*-model where they are 
distributed over higher energy states. Formulas includ- 
ing finite temperature corrections are derived in the Ap- 
pendix. 

In the limit of weak scattering, i.e. r+ -^ (tempera- 
ture dominated regime), self consistency is not required 
in Eq. I^J where Q on the right hand side can simply 
be replaced by lo. This corresponds to the temperature 
dominated regime with 7 <C T. Here 7 is the impurity 



scattering rate in the superconducting state and is fer- 
quency dependent. Considerable simplifications result as 
seen in the work of Hirschfeld et alm^ Generalizing their 
result for < 7 <C T to the non-equilibrium case, the 
real part of the conductivity is given by 



n^ r 
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where the scattering rate t^^{ijj) = 7(w) = — 5in(I;(a; + 
10+) is to be determined by Eq. Q with Cj replaced by 
w on the right hand side. The result for a general value 
of c is 



-^{uj)^ttT+ 



c^ + A+{uj) 



Aoc4 + 2c2^_(w)+v42_(w)' 



where 




A±{uj) = 


( ^"^ V 

VTrAo/ 



(f 



±ln 



/2Ao 



(9) 



(10) 



Two limits are normally considered. The Born limit with 
c ^ 00 and the unitary limit with c = 0. In these limits 
we get: 
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c — > 00. 



(11) 

In this last equation we can write 7rr+ = n//[7rA^(0)] 
and making use of the relation c~^ — 2TTN{Q)Vimp, we 
get 7rr+/c2 = 2rB with Tg = 2TTniN{0)Vg^p, the weU 
known formula for impurity scattering (Fermi Golden 
Rule). We can now use the scattering rates in the con- 
ductivity formula ||SJ). 

We begin with the /i*-model. For small temperatures 
<M*, (7<r</i*) 



ai{T,uj) 
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The finite frequency conductivity samples only the scat- 
tering time at the single frequency u = ^l* in this model. 
They are 



-1(m* 



7r-'r+ 



2r? 



2n (7r/2)2+ln2(^ 

2n, 







/2/n) ' 

c — > 00. 

(13) 
We see that in the Born limit r ^ (/x* ) is directly propor- 
tional to the square root of the excess quasiparticle den- 
sity, while it varies as its inverse in the unitary limit with 
some other additional, weaker logarithmic dependence on 
the square root of n. This means that the half width of 
the microwave conductivity will be strongly affected by 
changes in n and that the effects will be opposite in Born 



and unitary limit. These features are directly related to 
the details of the quasiparticle density of states at small 
Lo va. Si d-wave superconductor and t^^{ii*) as a function 
of /i* gives this information directly. 

From Eq. (|12|l one can easily compute the spectral op- 
tical weight remaining under ui. We define A as 
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duj ai{uj) 
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2n, 



(14) 



which gives a square root of n law independent of Born 
or unitary limit and of temperature, but the restriction 
7 ^ T ^ /i* must be noted. This accounts exactly 
for the missing superfluid density brought about by the 
excess quasiparticles. 

In the top frame of Fig. |31 we compare results based 
on the approximate Eq. (|12|l with full numerical results 
based on the exact expression ^. The same impu- 
rity parameters as presented in Fig. [3 are used, namely 
i+ = 0.1 meV with a gap Aq = 24 meV and T = 3 K. Sev- 
eral values of the non-equilibrium chemical potential ^* 
are considered: /i* = (black), ^* = 1.47 meV (dashed), 
^* == 2.94mcV (dotted), //* = 4.41 meV (dash-dotted), 
and ^* — 5.88 meV (dash-double dotted). The black 
curves are results based on Eq. (0) while the open sym- 
bols are based on Eq. H12|l . The frequency lu is restricted 
to small values below 1.0 meV. We see remarkable agree- 
ment between the full theory and Eq. H12|) . This agree- 
ment, of course, will fail as uj is increased outside the 
validity of the approximation used to derive Eq. l(T^. 
This is seen in the bottom frame of Fig. |31 where oj now 
goes up to 50 meV, i.e. up to about twice the gap value. 
For this higher frequency region no simple analytic re- 
sult can be obtained and it is necessary to return to the 
numerical results of the previous section. 

So far, except for Fig. [SJ we have emphasized the ex- 
pected frequency dependence of the optical parameters 
for a given value of excess quasiparticle density n. The 
data in Fig. [S] can be replotted to give the conductivity 
ai{uj^ fi*) at a given frequency lu as a function of fj,* . It is 
clear from Eqs. (|12|l and (|13|l that in the Born limit for 
ujT(fi) ^ 1, ai(T, uj) will be proportional to (/i*)^ and in- 
versely proportional to w^. This is shown in the top frame 
of Fig. where exact numerical results (solid lines) are 
compared with the approximate results of Eq. (|12|l (solid 
symbols). The expected linear variation with the excess 
quasiparticle density [n ex (m*)"] is pretty well verified. 
Since in experiments it is the change in reflectance R 
which is measured directly we show in the bottom frame 
of FiglHI results of this quantity at five frequencies, again 
as a function of n ex (/i*)^. These data, based on ex- 
act numerical calculations, show that the dependence of 
5R/ R on n is not exactly linear. 

Similar results to those embodied in Eqs. © to l(TH) 
for the /z*-model can be obtained in the T* model as 
well. We give only a few analytic results here without 
comparison to full numerical work. First we note that 
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FIG. 5: Comparison of exact results (black lines) for ai{T,u}) 
[Eq. ||2J] with the approximate analytic result of Eq. (I12II 
(open symbols). The top frame is for frequencies restricted 
to < w < 1 meV while the bottom frame spans to 50meV. 
The gray dashed curve is the normal state for f^ =0.1 meV 
and r = 3 K. The superconducting gap is 24meV, five values 
of /i* are shown. 



to 



0.5 



0.4 



0.3 



0.2 



0.1 



0.0 
3.0x10' 
2.5x10"' 
2.0x10' 
1.5x10' 
1.0x10' 
5.0x10' 

0.0 







co = 0.5meV ^^* 


' CO = 1 meV ^y^ 


■ , A , T ,Eq. (13) y^ 


/■ _ - ' "^ 
\/ ,-^-'' ^ 








(0 = 0.5 meV 


ffl = 1 meV . > '' 


co = 5meV •'■''' ^^^^^'^ 

(o=10meV ''^-'-^^ 

-■'' '^^-' 

■'^■■^t" 



10 15 20 25 

{^f (meV^) 



30 



35 



FIG. 6: Top frame: Real part of the conductivity (T\(y), n*) at 
three frequencies uj = 0.25, 0.5, and 1.0 me V as a function of 
the square of the chemical potential /i* (~ n). The lines are 
based on exact numerical evaluations of Eq. Q while the solid 
symbols are based on the approximate analytic formula 1121 . 
Bottom frame: The normalized change in refiectivity SR/R 
at five frequencies, uj — 0.25, 0.5, 1.0, 5.0, and 10.0 meV as 
a function of (/i*)^ (~ n) based on numerical evaluation of 
Eq. p. In both frames Aq = 24meV, t+ = 0.1 meV, and 

r = 3K. 



the residual absorption is 
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which is just the missing superfluid density due to T* 
[see Eq. Q]. The situation is, however, more compli- 
cated when the details of the microwave conductivity are 
considered. A simple form such as Eq. H12I) can no longer 
be derived and we must return to Eq. (jH)) which gives 
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where the temperature in the Fermi function is T* rather 
than T. Even if the frequency u in the logarithm that 
occurs in the Eqs. ^ for t{u) is replaced by T* , the 
integral in Eq. IjlGI) cannot be done analytically for a 
finite value of uj. In effect several values of t{uj) are 
sampled in Eq. ()16|l in contrast to the /i*-model where 
only r(/i*) enters. This is a true difference between the 
two models. If w = 0, however, we do get an analytic 
result (assuming 7 < T < T*): 



a{uj -^ 0) 



K 1 

47r 2Vb ' 



^i 



16n 



47r Tr'->T+ 



In^ 



c = 0. 



(17a) 



In the opposite limit, uj going to cx) 



constant value of 7 and find 
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We note that, in both cases, the dependence of the con- 
ductivity on the excess quasiparticle density n is very 
different in the two impurity limits considered. The na- 
ture of the impurity scattering is having a profound effect 
on the optical conductivity. For w — > in the Born limit, 
n drops out of the expression for the conductivity which 
then takes on its normal state value. For the unitary 
limit, instead, it goes like nln {n). For lj — > cxd in the 
Born limit the conductivity is proportional to n and in 
the unitary limit it varies inversely as In (n). Similar re- 
sults hold for the /z*-model; only the numerical factors 
are changed as can easily be verified from the general case 
(|12|l . In particular, the c ^ 00 result of H17a|) stays un- 
changed and the c — Q result is reduced by a factor of two, 
and the numerical factor in the logarithm tt/vo = 1.81 
is reduced by \/2 = 1.41; the w — > 00 result of Ijl7b|) 
is reduced by two and for c = we get the same result 
except for the numerical factor in the logarithm which 
goes from 1.81 to 1.41. While there are some quantita- 
tive differences between the results for the real part of the 
conductivity between T* and /x*-models, in the restricted 
parameter range in which our analytic results apply, and 
for w — > 00 the excess quasiparticles dominate the value 
of the conductivity. The condition 7 ^ T <C T* implies 
that n2/(4;r)(T*/Ao)2(2/3)7r2(rB/cj2) [which is the re- 
sult for the Born limit in Eq. Ijl7b|) ] is much larger than 
its value when 7i — Q, in which case T replaces T* in this 
expression; we note that the condition T <^ T* has been 
assumed. The same holds for the unitary limit. 

Finally, we consider possible modifications of the so- 
called universal limit. Applying a nodal approximation 
to the formula for the T < Aq conductivity, Eq. (O, it 
reduces at zero frequencji^ 



(20) 

The connection to croo depends linearly on n and on the 
inverse square of 7. Universality is lost. 



V. CONCLUSION 

We have presented results for the modifications of the 
optical conductivity ai{T,uj) that are brought about by 
introduction of a finite non-equilibrium excess quasipar- 
ticle density n. To describe the non-equilibrium dis- 
tribution in energy, we use two simplified models that 
have been found to be useful for the s-wave case. While 
not expected to be accurate the /x*-model of Owen and 
Scalapino2i and the T*-model of Parker^^i^^ have the 
great advantage that the physics involved becomes trans- 
parent. The frequency dependent conductivity shows 
distinct features associated with the chemical potential 
which should be observable. While both, the fx* and 
the T*-model give similar results qualitatively, there 
are some important quantitative differences. Numeri- 
cal results for the real and the imaginary part of a are 
given separately, as is the refiectance and the normalized 
change in refiectance brought about by the excess quasi- 
particle density. We have also given several analytic re- 
sults which we hope will prove helpful in the analysis of 
experimental data. 
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for T ^ this gives ai(0, 0) = [nll{^TT)][ll{TTl\o)] = aoo 
which is the well known universal limi t'^^i'^^ independent 
of impurity scattering. Generalization to non-equilibrium 
in the /i*-niodel gives 
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where 7(w) ~ 7[1 -I- &(a;/7)^] (with b a constan-feSi) ob- 
tained from a self consistent solution of Eq. Q for small 
but finite uj. For /x* <C 7 we can replace j{fJ.*) by its 



APPENDIX A: FINITE TEMPERATURE 
CORRECTIONS 

We consider finite temperatures in the penetration 
depth. The algebra is simplest in the T*-model. The 
thermal quasiparticle density in a d-wave superconduc- 
tor at temperature T is given by (7r^/12)(T/Ao)^ and so 
the excess non-equilibrium density n in the r*-model is 



TT 
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Thus, we have 
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which reduces to the known result T*/Ao = ^^firij-K for 
T = 0. 

It is easy to work out the smah T correction 
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while the small n correction to finite T is 
T* _ T 6nAo 



(A3) 



(A4) 



Equation ijSSll is vahd for Un/ir'^ > (T/Ao)^ and (|X4JI 
for 1211/71^ < (T/Ao)^. In (jA3|) the correction goes like 
T^ and is inversely dependent on y/n while in HA4(I it is 
linear in n and inversely dependent on temperature. The 
reduction in gap value due to the non-equilibrium density 
n also has a temperature dependence which is 
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which follows directly from Eqs. IJA3|I and ljA4p . We see 
from Eq. (|A7I) that the low temperature behavior of the 
penetration depth has changed from a T to a T^ law 
through the presence of the non-equilibrium density n. 
This is analogous to the effect of impurities which also 
bring about a T to T^ law change at low T. Of course, 
at higher values of T we recover a linear in T law as 
is indicated in Eq. IjASp with a small correction for the 
non-equilibrium density n. 

A similar situation holds in the fj,*-model but the al- 
gebra is not as tidy and we will only give two results. It 
is easy to work out the finite temperature correction to 
the chemical potential. From its definition (|A1I) 
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The solution for /i*/Ao is 



r21n(2). 
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which reduces to the known result ~32y/(3n)^/n^ at 
T = 0. It is easy to work out the lowest order T correc- 
tion to this result as well as at finite T, the lowest order 
n correction. As these corrections come in as higher or- 
der effects in the penetration depth we will not give the 
expressions here. 

The penetration depth follows directly from Eq. {T)) 
with T* replaced by Eq. HA2p . This gives 
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For T/Aq < \/V2JnpK^ (small T correction) we get the 
first correction for temperature 
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and for T/Ao > ^Ylnjir^ (small n correction) the first 
correction for non-equilibrium to temperature T is 



Ai (T*,T) Ai(0) 



l-21n(2)^-i^n^° 
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which properly reduces to \/2n at T = and to zero for 
n = 0. For ji* /T — > 00 the expression for A^ remains 
that given Eq. 10, and we get 
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and the square bracket reduces to 1 — -\/2n at T = and 
shows a positive leading, linear in temperature correction 
with the slope half the value it would have in the equi- 
librium case. This is different from the T*-model and is 
due to the very different non-equilibrium thermal distri- 
bution employed in the models which shows up at very 
small T. At large T, the /x*-model gives 
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Ai(r*,r) 
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l-21n(2)|- + ^Ao 
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which is the same as Eq. (IA8|) except for a slightly differ- 
ent numerical coefhcient in the last term 0.5 rather than 
0.8. 
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